Evidence of Raleigh-Hertz surface waves and shear stiffness anomaly in granular media 



L. Bonneau, B. Andreotti, and E. Clement 
PMMH, ESP CI, CNRS (UMR 7636) and Univ. Pans 6 & Pans 7, 
10 rue Vauquelin, 75005 Paris France. 
(Dated: March 27, 2008) 

Due to the non-linearity of Hertzian contacts, the speed of sound in granular matter increases 
with pressure. Under gravity, the non-linear elastic description predicts that acoustic propagation is 
only possible through surface modes, called Rayleigh-Hertz modes and guided by the index gradient. 
Here we directly evidence these modes in a controlled laboratory experiment and use them to probe 
the elastic properties of a granular packing under vanishing confining pressure. The shape and the 
dispersion relation of both transverse and sagittal modes are compared to the prediction of non- 
linear elasticity that includes finite size effects. This allows to test the existence of a shear stiffness 
anomaly close to the jamming transition. 

PACS numbers: 45.70.Ht, 05.45.Xt, 43.75.+a, 91.60.Lj 



Many disordered condensed systems - mostly multi- 
stable systems such as regular fluids trapped in a glassy 
phase PQ - can loose their shear rigidity, characteristic 
of ordinary solids, when they are on the verge to yield. 
For static, non-cohesive granular assemblies, the external 
pressure p is the only source of confinement that may jam 
the packing in the rigid phase. Thus, at zero thermal ag- 
itation, a solid/fluid transition would be reached under 
a vanishing confining pressure p i.e. when the packing 
becomes a marginal solid with just the minimal amount 
of contacts Z per grain suited to sustain a large scale 
elastic network, i.e. at isostaticity (Z = Z lSQ ). In the 
simple case of spherical, friction-less, disordered granular 
packing [2 , it has been shown theoretically that the local 
linear response becomes non-affine near jamming, due to 
the overwhelming appearance of weak floppy modes of di- 
vergent spatial extension [3J . These soft modes make the 
shear to bulk modulus ratio vanish at the transition. Can 
this anomalous linear response be evidenced experimen- 
tally? Is this jamming scenario robust enough to describe 
real granular matter, including frictional contacts [H[S], 
gravity loading or non-spherical particles, close to jam- 
ming? Does it provide a generic picture for the glassy 
transition of weak solids viewed from the jammed phase? 

In this letter, we show that acoustic surface waves 
sounding provides new insight into the structure of the 
elastic networks in the vicinity of jamming. This ex- 
perimental technique allows to measure accurately the 
elastic properties of a granular packing under vanishing 
pressure - namely, close to a free surface, under gravity 
loading. In ordinary elastic solids, surface waves (called 
Rayleigh waves) are a combination of compression and 
shear waves and travel at a speed slightly smaller than the 
bulk shear wave [5] (~ 5000 m/s for glass). In the context 
of granular matter, ethologists have reported that numer- 
ous species living at the desert surface use sound waves to 
probe their environement. Indeed, the first experiments 
of surface mode propagation were conducted by biologists 
who interpreted their signal as a collection of fast (bulk) 
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FIG. 1: (a) Experimental set-up. Acoustic source of sagital 
waves, (b) Acoustic source of transverse waves, (c) Signal of 
the accelerometer mounted on the source, (d-e) Transverse 
wave-packets received at x = 60 cm from the source, for two 
preparations of the granular packing (see text). Note the shift 
of the wave-packet centre (•) and of the phase o. 



waves and slow (shear) waves although the measured 
speeds (~ 50 m/s) were incredibly small j7]. This was 
confirmed by field measurements performed at the sur- 
face of sand dunes in the context of booming avalanches 
studies [SI IS] ■ Recently, these surface waves were ascribed 
to the gravity induced index gradient, which plays the 
role of a wave-guide, and analysed in two slightly differ- 
ent models of non-linear elasticity [TOj [IT] . 

To clarify the theoretical issues addressed here, let us 
rephrase in the framework of non-linear elasticity the pre- 
dictions of Wyart et al. [3J on the elastic anomaly in- 
duced by soft-modes. For geometrical reasons, the Hertz 
contact force between two grains depends non- linearly 
on their relative inter-penetration 5. On this basis, the 
macroscopic elastic free energy \§\ of an isotropic granu- 
lar packing can be written in the very general form |14j : 



T = E 



J :BS 5 / 2 



(1) 



2 



is the coarse-grained strain tensor, 8 



where 

— Tr(ujj) is the volumic compression and u®j = tiy + gJij 
the traceless strain tensor. A and B are two dimension- 
less elastic coefficients that characterise the material stiff- 
ness under shear and compression, respectively. We as- 
sume that the average number of contacts per grain Z is 
sufficient to characterise the microscopic packing geome- 
try and thus, that A and B are function of Z. Indeed, in a 
frictional packing, different values of Z can be obtained 
under the same pressure p; Z and p are thus indepen- 
dent state variables [5]. It should be emphasised that 
T is not supposed to describe the stress-strain curve ob- 
tained from a loading test, which is composed by series 
of elastic loadings at fixed Z and of plastic events during 
which Z changes. 

The mean field calculation [11 predicts that A and B 
vary linearly in Z and thus remain finite at p = 0. On the 
other hand, the jamming theory predicts that upon ap- 
proaching isostaticity, the stiffness ratio A/B exhibits an 
anomalous behaviour - as well as many other quantities - 
and scales on the excess of contacts above the isostatic 
value, Z — Z- lso cx S 1 ^ 2 . So, this approach predicts that 
A/B vanishes with pressure. Let us emphasise again that 
the reality of the jamming point as a critical point was 
assessed for frictionless soft spheres but there are several 
indications that this line of ideas could be generalised to 
frictional packings [HE3 an d even to glasses [TB] , 

Starting from expression ( II) , the surface elastic modes 
can be derived as shown in [XX] . They propagate through 
an infinite though discrete collection of modes polarized 
sagitally or transversally. The mode labelled n of wave- 
length A penetrates the sample over a typical depth nX 
and thus feels a typical pressure p = pgnX. Then, the 
mean field theory predicts a velocity dependence of the 
form: v <~ (nA) 1 / 6 while the "soft modes" theory pre- 
dicts, instead, v ~ (n A) 1 / 3 . Thus, the dispersion relation 
of surface waves constitutes a direct experimental test of 
the existence of a shear stiffness anomalous scaling. 

As confirmed by preliminary experiments [5], the weak 
dependence of the speed of sound on n makes the ex- 
periment very difficult to control and almost impossible 
to analyse. Indeed, using any standard source, a huge 
number of propagating modes are excited that remain 
superimposed over large distances. To bypass this prob- 
lem, we have successfully designed an experiment that 
isolates the first sagittal and transverse modes. For this, 
measurements are performed in a rectangular channel of 
width W — 20 cm and length 180 cm, which serves as a 
second wave-guide (Fig.JTJi): the lateral boundary condi- 
tions impose a relation between the measured wavelength 
A and the wavelength Aoo that would be selected in an 
infinitely large channel, at the same frequency: 
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FIG. 2: Space time diagram showing the wave-packet prop- 
agation. Starting from the raw signals received at different 
positions, the wave-packet is roughly localized by computing 
the signal envelope. Then, a local fit by a gaussian wave- 
packet allows to determine the centre of the wave-packet (•) 
and its phase with respect to the source (symbols ■ show the 
space-time coordinates of an iso-phase event). The best fit 
(thin lines) allows to extract the group and phase velocities. 



Moreover, the sources are conceived and tuned to ex- 
cite essentially the modes n = 1: sagittal waves are pro- 
duced by an electromagnetic shaker (without any spring) 
whose axis is finely guided by a ball bearing slider cou- 
pled to a very rigid transverse metallic blade; transverse 
waves are produced by a rough cylinder inside which a 
permanent magnet vibrates under the action of a mag- 
netic field (fig. [ljj) . The channel is filled of glass beads 
(E = 70 GPa, d = 150 //m) over the height H = 20 cm. 
For such an aspect ratio, with smooth boundaries, the 
Janssen effect is negligible so that the pressure is ex- 
pected to vary linearly in depth. The acoustic isolation 
is insured by 20 cm thick boundaries. 

The experiment is also designed to prevent another 
problem. As the sample presents random heterogeneities, 
the acoustic signal is composed by an effective medium 
response and a coda related to speckle effect [12]. Their 
relative amplitude is controlled by the number of grains 
in contact with the transducer. We have chosen to work 
with accelerometers of diameter D = 13 mm which al- 
low to measure the three components of the accelera- 
tion in the bulk of the sample. Around 3 10 4 grains are 
in contact with the transducer and the measured ampli- 
tude of the coda-tail is around 5% of the coherent signal. 
By comparison, other techniques like a laser vibrometer 
would only probe the rough surface of the packing and 
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FIG. 3: First sagittal mode, (a) Amplitude and (b) phase of 
the vertical (•) and horizontal (o) displacement as a function 
of depth, for / = 315 Hz (A ~ 21.5 cm), (c) Length A over 
which the vibration decays as a function of the wavelength A. 

due to the small size of the spot (5 mm) would only av- 
erage over ~ 400 grains, yielding a coda and a coherent 
signal of the same order of magnitude. We have also 
checked that the propagation is not affected by the pres- 
ence of other accelerometers between the source and the 
receiver. Besides, the transducer size D should be at 
least a fraction of the wavelength A, which imposes to 
work at rather low frequencies / < 1 kHz. In summary, 
the experiment has to be analysed having in mind the 
hierarchy of length scales: d <C D < A < H. 

The typical vibration amplitude we use is ~ 10 nm, 
although the propagation properties remain the same up 
to ~ 100 nm i.e. a strain of 10 -6 . Above, new peaks ap- 
pear in the signal and period doubling is observed in the 
coda. The preparation of the sample is amongst the most 
difficult part of the experiment. Prior to each measure- 
ment , we systematically sweep a thin blade longitudinally 
and transversally through the packing in order to remove 
any memory effect due to the granular initial filling or the 
subsequent accelerometers manipulation. The extra-sand 
above the level of the channel is then gently removed, let- 
ting a flat surface. The figure [l] compares a typical signal 
obtained with this procedure (d) to that obtained when 
pouring the grains and flattening the surface by tapping 
with a hammer on a plaster float (e). In the first case, 
the phase and the travel-time of the wave-packet is re- 
producible. By contrast, the apparent phase velocity and 
group velocity can vary up to 25% from one sample to 
the other, although they are macroscopically identical. 
This non-universality is the first important conclusion of 
this letter: the packing elastic properties depend on the 
the preparation protocol. 

Fig. [2] shows the positions of the centre of the wave- 



packet and of an iso-phase event as a function of time. 
For this, we prepared ~ 20 realisations of the packing. 
For each, the signals of four accelerometers placed at 
different positions are acquired and fitted by a gaussian 
wave-packet. The propagation of a coherent mode - and 
not speckle - is clearly evidenced by the linear relation- 
ship between space and time and by the reproducibility 
over independent microscopic realisations. The slopes of 
the relations respectively give the group and the phase 
velocities, v g and v v . Looking now in depth, the vibra- 
tion amplitude is observed to decrease over a distance 
A, of the order of a half wavelength (fig. |3j;). The sagit- 
tal waves are elliptically polarized, with their principal 
axis in the vertical direction and along the direction of 
propagation. The figure [3] compares the shape of the first 
mode to the prediction of our model [TTj . We obtain a 
fair agreement with the theory which confirms that the 
first mode has been successfully isolated, as desired. 

In our set-up, the working frequency / can be varied in 
the range 300—750 Hz for sagittal waves, and 300—550 Hz 
for transverse waves. Outside these frequency ranges, the 
signals received are either distorted due to the source, 
or multi-modal. No propagation at all is observed below 
200 Hz. We have measured v g and v v , as well as their sta- 
tistical uncertainty, every 3 Hz. These values are used to 
reconstruct accurately a single dispersion relation /(A -1 ) 
(fig. [4]) that simultaneously fits the group and phase ve- 
locities in the least square sense. As expected, the prop- 
agation is dispersive since v g and v v are different. This 
can be related to the two wave-guide effects previously 
mentioned. Due to the finite depth and width, a cut- 
off frequency appears in the dispersion relation, which 
corresponds to the first resonance of the system. The 
group velocity is expected to vanish at this frequency 
and the phase velocity to diverge. This explains the be- 
haviour observed at low frequency, and in particular the 
increase with / of the ratio v g /v v . At high frequency, 
an asymptotic behaviour controlled by the pressure in- 
duced wave-guide is reached in the limit of wavelengths 
A small in front of the channel transverse dimensions (H 
and W). The ratio Wg/w^ tends to a constant equal to 
0.82 ± 0.04 for the preparation described above. As v g is 
the slope of the dispersion relation, v g /v v is the scaling 
exponent between / and A -1 . The measured value is very 
close to that expected if A does not vanish at the surface 
(5/6). Thus, our experimental results fully confirm the 
Hertzian picture down to A ~ 250 d and does not show 
any evidence of anomalous exponent (2/3) when p — ► 0, 
that would be associated to the jamming point [3J ITS] . 
It would be interesting to pursue this technique to lower 
the probing wavelengths in order to push the limits to 
weaker confining pressures. But to this purpose, the lim- 
itations due to finite probe size, preparation sensitivity 
and speckle noise have to be overcome. 

Globally, one can compare the experimental data to 
the prediction of the model, assuming that A does not 
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FIG. 4: (a) Dispersion relation of sagittal (o) and longitudinal waves (A): frequency / as a function of the rescaled wavenumber 
H/\. (b) Corresponding group v g (white symbols) and phase (black symbols) velocities, (c) Ratio of the group and phase 
velocities v g /v v as a function of H/X. In the three graphs, only one fifth of the measured points are shown. The predictions of 
the model, including the finite width and depth, are shown in solid (sagittal) and dotted (transverse) lines. 



vanish at the surface and remains nearly constant at the 
scale A and taking into account the finite width and depth 
(lines on fig.|4|. The agreement is excellent (within 5%). 
The dispersion relations of transverse and sagittal modes 
turn out to be nearly equal (fig. |4J|. This striking be- 
haviour is one of the robust output of our model 
For ratios B/A = 0(1), the difference would only be of 
10% and would be almost indistinguishable when, say, 
A < 0.2 B. The physical reason is that the restoring 
force for both modes is the shear elasticity (parameter 
A). In the limit A <§C H the dispersion relation of the 
first modes takes the form: 

/ = a(£;/p) 1/3 5 1/6 A- 5/6 with a ~ 0.77 A 1/2 B~ 1/6 (3) 

The best fit gives: /A^ 6 — 77 ± 1 m 5 / 6 s _1 , which cor- 
responds to a value of A 1 I' 1 B~ X ^ ~ 0.23. By contrast, 
the mean-field expectation is 0.40 for frictionless grains 
and 0.61 for infinite friction. Thus, the measured shear 
stiffness is 3 to 5 times smaller than expected from the 
mean field theory, as observed in numerical simulations 
|13j . The "soft-modes" theory is the only one explaining 
this mean-field failure. Still, to be consistent with our re- 
sults, one has to conclude that the packing does not tend 
to isostaticity in surface: for frictional packings prepared 
in a simple way, Z remains significantly larger than Z mo 
under vanishing pressure. 
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